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Abstract 



We give the details of symplectic quantization for a system containing second class 
constraints. This method is appropriate for imposing infinite series of constraints 
due to the boundary conditions. We use this method for massive bosonic strings in 
a background B-field and find the correct expansions of the fields in terms of the 
O ' physical modes. We have found a canonical basis for this model. 
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1 Introduction 

oo 

Quantization process has been one of the most important problems from the early days of 
quantum theory. According to the famous prescription of Dirac [1], given a classical theory 
with a well behaved bracket in the phase space, the corresponding quantum theory comes 
out as the result of changing the dynamical variables to operators. Hence, the consistent 
algebraic structure of classical brackets is an important part of the quantization process. If 
the system possesses second class constraints one needs to consider Dirac brackets instead 
^ ; of Poisson brackets. However, this procedure should be consistent with the dynamics of the 

system both in the classical and quantum levels. This demand makes us to investigate the 
consistency of constraints in the course of the time and finally find a consistent algebra of 
weakly vanishing brackets among the constraints and the Hamiltonian[2l [3]. 

On the other hand, we have difficulty with the Dirac method since by restricting the 
system to live on the constraint surface one may miss the beautiful Poisson structure of the 
phase space. In other words, we need a new canonical structure that guarantees a closed 
algebra of Poisson brackets on the constraint surface. According to the famous "Darboux 
theorem" one can, in principal, find a set of canonical coordinates in which the constraint 
surface is described by a number of canonical conjugate pairs in addition to some extra 
coordinates. Faddeev and Jackiw jl] used the above theorem, at each step of consistency, 
to construct a new method for analyzing the constrained systems. The Faddeev- Jackiw 
formulation and its equivalence with the Dirac approach is studied in detail by several 
authors, see for instance |6]. 
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A modified version of Faddeev- Jackiw formulation recognized as " symplectic quantiza- 
tion" is more or less used by some authors, when they want to build a quantum theory 
out of a given classical actionjTJ El O EE]- As we will see, the most essential point in this 
process is finding the complete set of physical modes and the correct expansions of the fields 
in terms of them. The main goal of this paper is giving the essential aspects of symplectic 
quantization and applying it for quantization of massive bosonic string in a background 
B-field as a rich example. We will show that the method in its original form works well 
for the bosonic string in background B-field with no need to any modification such as the 
additional time averaging suggested in [9]. 

In this paper we first review the essential aspects of the symplectic quantization in 
section 2. We restrict ourselves to a purely second class system. This is specially the case 
when all the constraints have originated from the boundary conditions. In section 3 we 
consider in details the massive bosonic string again in a background B-field. As we will see 
the constraint structure of this case is more than a simple extension of the massless case. 
Then we exhibit the power of the symplectic quantization method in this interesting and 
illuminating example. The results on the commutation relations would be discussed in the 
text and also in the last section which is devoted to the concluding remarks. 



2 The Symplectic Approach 

The Hamiltonian formulation of a dynamical system can be achieved by using the first 
order Lagrangian 

L F .o. = ^PiQt ~ H(q,p). (1) 

i 

If one uses general coordinates (i = 1, 2N) for the phase space, which may be or may 
not be canonical, the first order Lagrangian would be written as 

Lf-.o. = 5>(0f -#(£)■ (2) 

i 

In this form the action is 

S = J [a^)de - H(£)dt], (3) 

where the kinetic term exhibit a one-form on the phase space manifold. The Euler-Lagrange 
equations of motion can be written as = dH/d^ 1 , where Uij = dat/d^ — dcij/dl; 1 is 
the symplectic matrix. It defines the "symplectic two- form" as 

Q = -Uijdg A dg = da, (4) 

where a is the kinetic one-form of the action ([3]). The equations of motion can be solved 
as £ l = u^dH/dtf , where uj^ is the inverse of Uy. These equations should be equivalent 
to the canonical equations of motion £ l = H}, where { , } means the Poisson bracket. 
Hence, for any two functions /(£) and g(£) the Poisson bracket should be defined as 
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"ab = — TZT^r (6) 



As we see, the important role of the symplectic two-form is that the inverse of symplectic 
matrix defines the Poisson brackets in arbitrary coordinates of phase space. 

Now, suppose we are given a set of second class constraints (p a (0 ~ on the phase space 
variables. Suppose this set includes primary constraints as well as secondary constraints 
which come out from the consistency of primary ones. In other words, no new constraint 
may be added to the system as the result of consistency of the present constraints. For 
instance, these constraints can be considered as the infinite set of boundary conditions 
and their consistency conditions, to be discussed in the following section. Under imposing 
the constraints on arbitrary coordinates £ fc of the original phase space, the reduced phase 
space is described by the coordinates rf a — 1, .., 2m (m < N), which are not necessarily 
canonical. On the reduced phase space, we have £ fc = £ fc (?7 a ), and the induced symplectic 
tensor is 

drf drf 

The induced symplectic tensor u a b can be written by imposing the constraints on the 
symplectic two-form Q. Inverting cj a b, then gives the inverse tensor u ab which determines 
the Dirac brackets on the reduced phase space. 

For a field theory in d + 1 dimensions with the dynamical fields (j) s (x,t) the symplectic 
two-form in terms of the original fields is given by 

n = / ddxdu s(x, t) A d(f) s (x, t). (7) 

s J 

where H(x, t) are momentum fields. Upon imposing a set of second class constraints, such 
as initial boundary conditions and their consistency conditions, suppose the phase space 
fields 4> s (x, t) and Tl s (x, t) can be written in terms of a restricted set of physical mods a n (t). 
Inserting the corresponding expansions of 4> s (x,t) and U s (x,t) in Eq. ([7|), the symplectic 
two-form can be written, in principal, as 

VI = J2 UJ nmda n Ada" 1 . (8) 

n,m 

The summations in (jSJ) are understood to include integration over the continues variables 
in cases where physical modes a n (t) depend on continues parameters. The inverse tensor 
u™ defines the Dirac brackets of the reduced phase space coordinates as 

{a n , a m } = u nm . (9) 

Note that the initial Poisson brackets is no longer valid on the reduced phase space. 

We emphasize the important point that up to this point the symplectic quantization is 
done without solving the equations of motion of the variables or fields. In other words, we 
have used the first part of the action ([3]) so far. The dynamics is, however, the responsibility 
of the second part, i.e. the Hamiltonian. Properties of the Poisson structure on the classical 
phase space, as well as the algebra of the operators in quantum mechanics do not depend 
at all on the form of the Hamiltonian. For example for a particle in one dimension the 
most important thing for quantization is the classical bracket {x,p} = 1 or equivalently 
{a, a*} = j^; which converts to [x,p] = ih or [a,a^] = 1 upon quantization. These relations 
do not depend on the explicit time dependence of a(t) and a^(t) to be for example a(0)e tuJt 
and a^(0)e~ lult for simple harmonic oscillator. 
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It worths reminding that upon ordinary conditions, the classical brackets, as well as the 
quantum commutators, are independent of time. To see this explicitly, suppose at time t 
we have {£ l (t), ^{t)} = uj % K Assume u) % i is independent of £ l -s, which is the case in most 
of the problems. In the course of the time, the variation of £ l is S^ 1 = H}5t. Using the 
Jacobi identity the bracket of and £ J at time t + St is 

{e+se,?+5e} = ^ +st[{{e,H},e} + {?,{?, h}}] 

= u* + 6t{H,{?,e}} 

= (10) 



3 Massive bosonic string 

Consider the Lagrangian of the massive bosonic string in an external B-field, 



L — - i 
2 J 



i 

12 



X 2 - X' 2 - m 2 X 2 + 2B ij X i X! i 



where "dot" and "prime" represent differentiation with respect to r and a respectively and 
X 2 = XiXi, etc. This is the simplified version of a model in which among the whole set of 
bosonic fields X^ an even number X % are influenced by a background antisymmetric B-field 
and are attached at the end-points to D-brains. Here, we have omitted those components of 
X M which possess ordinary Neumann boundary conditions and are not coupled to B-field. 
We consider the simplest case where i — 1,2, and the metric of the truncated target space 
is Euclidian. Hence, we consider all the space indices as lower indices. The antisymmetric 
B-field in two space dimensions can be written as 

Mi!)- (12) 

The massless case (m = 0) is studied in different aspects by several authors [Hi [12j [13], [TJJ 
[T5] with the well-known result of non commutativity at the end-points. Demanding the 
variation of the action vanish under arbitrary variation SXi gives the equation of motion 
as (d 2 — d 2 — m 2 ) Xi = 0, while the boundary conditions are X/ + BijXj = at a = and 
a — I. The momentum fields are given by Pi = Xi + B^Xj. The canonical Hamiltonian 
reads 

i r l 

H=- da [(P - BX') 2 + X' 2 + m 2 X 2 ] . (13) 
The fundamental Poisson brackets (before imposing the constraints) read 

{X l (a,r),X J (a',r)} = 0, 

{P l (a,r),P J (a,T)} = 0, (14) 
{Xi(a,T),Pj(a,T)} = SijS(a - a'). 

Consider the boundary conditions as primary constraint $^ = 0^| (7= o and $^ = (pi \ a =i 
where 

4 1] = MijX'j + B l3 P 3 . (15) 
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in which M = 1 — B 2 . As indicated in details in references [TTJ [12], consistency of primary 
constraints in time determines the Lagrange multipliers \ and Aj in the total Hamiltonian 
Ht — H c + Aj<&£ + Ai < 5 l - 1 ' ) to be zero and at the same time gives second level constraints 



as $^ 2) = ^i(0,r) and = ipi(l,r) where 



(2) 



^i(er, r) = c^P; - m 2 B ij X j . 



(16) 



Direct calculation shows that third and forth level constraints are (d 2 — m 2 )0j and (d 2 — 
m 2 )ipi at the end-points, respectively In this way even and odd level constraints at a = 
read respectively as 



(2r) = 
i 

(2r+l) 



(di-m'Y'^i U=o~0 r>l , 
= (dl - m 2 y<\> % | CT=0 ^ r > , 



(17) 



with similar expressions for $f r) and $f r+1) at a 



I. 



Expanding the real fields X(a, r) and P(<j, t) as Fourier integrals, we have 

1 



Xi(<T,r) 



Pi(a, r) 



1 



2tt 



oo 
oo 



dk[ai(k, t) cos ka + bi(k, r) sin fccx], 
dA;[cj(fc, r) cos A;a + cZj(fc, r) sin fcer]. 



It is clear that a,i(k, r) and Ci(k, r) are even while bi(k, r) and r) are odd functions of k. 
Before imposing the constraints all of the above modes are present. Using the fundamental 
brackets (114j) one can find the following Poisson brackets among the Fourier modes 



{di(k, r), Cj(k', r)} = 8ij8(k - k'), 

{b l (k,T),d j (k',T)} = 5 lj 5(k-k'), 



(19) 



while all other bracket vanish. This shows that the canonical pairs (aj(fc, r), Ci(k, t)) and 
(bi(k , t) , di(k , t)) act as an alternative canonical basis for the original phase space. Upon 
imposing the constraints a large number of these modes will be omitted, leaving a much 
smaller number of them as the canonical coordinates of the reduced phase space. 
We now impose the set of constraints ( ITT)) on the Fourier expansions (ITS)) to find 



dk(-l) r (k 2 + m 2 ) r [BijCj(k, t) + kMijbjik, t)] = 



/oo 
dk(-l) r (k 2 + m 2 ) r [kdi(k, t) - m 2 B ijaj (k, t)] = . 
-oo 



(20) 



Since the expressions in the brackets are even functions of k, the conditions (I2TJ)) for arbitrary 
r will be satisfied only if b(k, r) = —M~ 1 Bc(k, r)jk and d(k, r) = m 2 Ba(k, r) / k. This leads 
to omitting the modes b(k, r) and d(k, r) in favor of a(k, r) and c(k, r) respectively. For 
simplicity we omit the indices • • • on the fields and modes from now on. Imposing the 
constraints at the end point a = I gives 



dk(-l) r (k 2 + m 2 ) r 
dk(-l) r (k 2 + m 2 ) r 



-m 2 B 2 



kM 



-k + -m 2 B 2 M- 1 
k 



a(k, t) sin kl = , 
c(k, t) sin kl = . 



(21) 
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These equations show that a(k,r) and c(k,r) should vanish except for k = nix /I with 
integer n or when 

m 2 B 2 l2 

-k 2 Q . (22) 



k 2 = m 2 B 2 M~ ] 



1 + B 2 

The first possibility leads to oscillatory modes as 



n=l 

oo 

£ 

n=l 



a (n) M cos — M _1 5c (n) M sin — — 



c (n) (r) cos 



/ nn 
nna I 



I mi 



-m 



nno 



r sin ■ 



/ 



(23) 



The normalization coefficients \j2/l are chosen for future convenience. This choice also 
makes the correct dimensionality for as (Mass)" 1 / 2 and for as (Mass) 1 / 2 . Note that 
the fields Xi are dimensionless. 

The possibility (122]) corresponds to zero mode solutions with sinh k Q a and cosh k^a. 
Traditionally the zero mode solution is denoted as the zero frequency (infinite wave length) 
limiting term in the Fourier expansions, as shown for example for the massless case of the 
current problem in [TT]. Here, however, we interpret the zero mode solution as a solution 
which satisfies the boundary conditions not only at the end-points but also throughout all the 
medium. This interpretation works well for constant term in the case of ordinary Neumann 
boundary condition of free bosonic string as well as for the zero mode terms for the massless 
string in B-field. To see the details, suppose the first and second level constraints are valid 
throughout all the string as the following coupled deferential equations 



BP + MX' = . 
P' - m 2 BX = 



(24) 



which gives MX" + m 2 B 2 X = . Since B 2 
solutions of Eqs. ff24|) can be chosen as 



-B 2 1 and M = (1 + B 2 )l the most general 



X zm (cr,T) = — 



1 

1 



a (0) (r) cosh[A; (a - 1/2)} - —M^Bc^ir) sinh[fc (a - 1/2)} 

ko 



P zm {o-,T) = — 



r) cosh[k (a - 1/2)} + ^m 2 5a (0) (r) sinh[A;o(<7 - 1/2)} 



(25) 



where ko is given in Eq.( )22|) . Again we have imposed the normalization coefficient a/1// 
for making correct dimensionality and future convenience. Note that we have chosen the 
argument of sinh and cosh functions as measured from the center of the string. Besides 
respecting the symmetry of the Hamiltonian under a —> —a when we shift the origin of a 
to 1/2, the reason for this choice is that the physical modes in this sector turn out finally 
to be canonical modes, as we will see in the near future. The most general solution of the 
fields are superposition of zero-mode and oscillatory solutions as 



X = X Z 
P = P, n 



- X 

P 

OS- 



(26) 
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Now we want to quantize the theory using the symplectic approach. After a long but 
direct calculation the symplectic two-form Q = £\ J dadPi(a, r) AcOQ(er, r) in terms of the 
physical modes and c'"'(t) read 



n=l 



„ = smhW ( icf) A ^ + £ (i + ) (<fc<»> A ) (27) 



The important point in derivation of Eq.( l27j) is that all cross terms composed of zero- 
modes wedge oscillatory modes have been canceled. This is a good news, since it makes the 
symplectic matrix block diagonal with one 4x4 block for the zero modes and an infinite 
dimensional block for the oscillatory modes. Moreover in oscillatory modes there is no 
mixing between differen oscillators (i.e. different n). In zero mode sector also we do not 
encounter cross terms such as da^ A da^ and dcf^ 1 A dcf\ This is because of our suitable 
choice of zero modes in Eqs (]25|) as combinations of sinh[/c (o" — l/ty] an d cosh[/c (cx — 1/2)]. 
These opportunities make it easy to find the inverse of symplectic matrix and write down 
the brackets of physical modes. Hence, the non-vanishing Dirac brackets among physical 
modes follows 

[a ( T\c ( j s) ] = N- 1 5 l3 5 ns , (28) 

where 

si nh(fc Q N _ 1 k 2 l 2 

k l n 2 7l 2 

Note that we did not need to solve the equations of motion up to this point. The theory can 
be quantized, without any need to full solutions of the equations of motion by converting 
the brackets (128]) into the quantum commutators. Similar to the massless case, one can find 
the time dependence of the physical modes by directly solving their equations of motion. 
For this aim it is more economical to write the canonical Hamiltonian f fT3|) in terms of the 
physical modes. The result is 



JV = =^, Ar B El + ^L^0. (29) 



CO 



i n=0 



M-^] 2 + Mulla^] 2 , (30) 



n 2 7r 2 



where 

ul = m 2 M , u 2 n = m 2 + n ^ (31) 

The Hamiltonian (|30|) is a superposition of infinite number of independent harmonic oscil- 
lators. The canonical equations of motion can be solved for n = 0, 1, • • • as 

a W( r ) = a in) (0) cosu n r + -^—c {n) (0) smu n r , 

c (n) ^ = c (n) (q) cog ^ _ MuJna (n) (q) gin ^ ? (33) 

where a (n) (0) and c (n) (0) are Schro dinger modes. Inserting: and c^(t) from Eqs. 

( )32|) in the expansions fl23|) and fl25j) of the fields, gives in fact the solutions of the equa- 
tions of motion consistent with the boundary conditions, written in terms of the canonical 
Schrodinger modes. 

Our results here are different from those of reference [10] in the following aspects: 
1) The algebra of physical modes and c- n > (r) is much simpler. They are canonical 

conjugate pairs with position-momentum like brackets given in (|28|) . 
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2) Inserting our final expansions of the fields in terms of the Schrodinger modes in the 
symplectic two-form leads again to expression (|27|) for the initial values of the modes. In 
other words, there is no explicit time dependence as expected. However, direct calculation 
shows that, similar to the massless case of reference [9], the symplectic two- form constructed 
from expansions given in [10] contains time dependent terms. 

It seems that real functions for spatial dependence of the fields should be accompanied 
naturally with real time dependent of the modes as seen in ( 1321) . Combining spatial real 
functions with the time dependence of the form e tulT may lead to unwanted time dependence 
in the symplectic two- form. In fact, this time dependence is the origin of an unnecessary 
step of time averaging of the symplectic two- form suggested in reference [9]. 

In order to complete our discussions and enlighten some unclear points of the literature 
let us calculate the Dirac brackets of the original fields to see the effect of the B-field and 
boundary conditions on the commutativity of fields. To do this using the brackets ( 1281) 
and Eqs.( l23p .( 125]) and ( 126|) we can compute the equal time Dirac brackets of the original 
coordinates and momentum fields as follows 

[X i {a,T),X j (a',T)] = 2(BM-%f(a + a') , (33) 
[Pi{a, r), P,V, r)] = 2m 2 B, l3 f(a + a') , (34) 
[X i {a,T),P j (o',T)] = 6 ij g(<T,o') (35) 



where 



sinh knl ir ' n \ n 2 Ti 2 J \ I 

n=l 



and 



/a o, cosh[fc (<r + a' - I)} 

g{a,a) = 2k — — 

sinn k l 



4 l 2 kn \ nucr mra' 



COS — : — COS 



I ^— ' V ' n 2 7T 2 
n=l 

4 ^ i 2 ^ ( \ _1 n7ra n7Ta f 



71=1 



-TT 1 + sin — sin— . 37 

. n 2 ir 2 V n 2 ir 2 J I I 



Let us first consider the brackets (|33|) and (|34"|) in details. The function f(a) reduces, for 
m — 0, to 



„ . . a 2 1 /"nira\ , , 

/oW = - 1 + _ + _^_ sln (_) , (38) 

n=l 

which gives 



{X l (a 1 r),X j (a',r)} = a,a'^0, 

{XiiO^Xj^r)} = -2{M- l B)^ (39) 
{X i (/,r),X i (/,r)} = 2(A-r 1 J B) iJ , 

which is the standard results of non commutativity of the end-points of a massless string 
in the background B-field. For m^fl using the Fourier expansion 

, r, / ,m 2sinh£;oZv^ 1 ( n 2 n 2 \ /nna\ , A . 

**M„ - 1)] = E - ( ^ripi J - (— ) • («) 

n=l 
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in the interval [0,2/], it is easily seen that f(a) vanishes for every value of a in the open 
interval (0,2/). For a — and a = 21, however, the Fourier expansion ( 14"0"|) is no longer 
valid since, beyond the interval (0,2/), the hyperbolic function on the left hand side of 
Eq. (T40|) can not be expanded in terms of periodic functions sin (nira/l). Similarly, the last 
expression in the right hand side of Eq. (138]) is the expansion of the non-continues sawtooth 
function with non-continuity on the points a = 2kl for integer k. At the end-points a = 
and a = 21 the function fo(a) is +1 and -1 respectively from its definition in Eq. ( 1361) . The 
above argument can be used exactly in the same way for f(a) given in Eq. ( 136]) . In other 
words, f(a) in Eq. ( 136]) . for arbitrary k , is the same as given in Eq. fl38|) for the case 
ko = 0. Hence we have 



f{a) = for < a < 21, f(Q) = -1, /(2/) = 1 



(41) 



This result shows that non-commutativity of the coordinate fields at the end points are 
exactly the same for massless and massive cases. 

As we see from Eqs. ( 133]) and (134"]) . for m ^ 0, the momentum fields as well as the 
coordinate fields are noncommutative at the end-points. This is due to dependence of mo- 
mentum fields to coordinate Fourier modes a\ in contrast to massless case which is just 
composed of momentum Fourier modes cf 1 . It worth remembering that mixing of coordi- 
nate and momentum Fourier modes in the expansions of X(a, r) and P(cr, r) is basically 
due to mixed boundary conditions which is in turn resulted from the presence of B-field. 

Next we consider the bracket of a coordinate and a momentum field in Eq. (|35|) . The 
function g(a, a') given in Eq. (|37j) can be written as 



4 
7 



cos 



where 



g(a) = 2k 



E 

n=l 

cosh[k (a — /)] 



nito nircr 

- COS ; 



g(a + a 



sinh kol 



E 



n=l 



l 2 k\ 



I ^-f \n 2 7T 2 + l 2 kl 



nno 



cos 



Using the Fourier expansion 



cosh[A;o(cr — /)] 



sinh kol 
knl 



Ei 2k§l sinh kol \ (nita 
I I COS 



n 2 7r 2 + kll 2 



n=l V l 

for cre[0, 21], it is easily seen that g(a) = 2/1. So we have finally 



V / 



(42) 



(43) 



(44) 



2 4^ 
g(a, a') = - + - ^ 



/ / 



nna rnra 
cos — - — COS 



n=l 



I 



I 



5 (a - a') 



(45) 



The interesting point is that the result is again independent of m or ko- In other words, the 
fields X(a, r) and P(a, r) are still canonical conjugate pairs in the reduced phase space, for 
the massive case. 



4 Concluding Remarks 

Our main goal in this paper is to reintroduce "the symplectic quantization method", in 
such a way that is appropriate for imposing the second class constraints originated from the 
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boundary conditions. We showed that the essential point is imposing the set of constraints 
on some "appropriate expansion" of the fields in terms of suitable modes. Although the 
full dynamics of the system is not required in order to use the machinery of the symplectic 
approach, it is necessary to investigate the dynamics of the constraints; which is, in fact, 
the program of finding the constraint structure of the system. 

The main outcome of this procedure is recognizing all independent physical degrees 
of freedom, i.e. the coordinates of the reduced phase space, and finding expansions of 
quantities of interest, such as the original fields and the Hamiltonian in terms of them. Note 
that in this viewpoint the mode expansions of the fields are not just combinations of the 
solutions of the equations of motion with some meaningless coefficients. These coefficients, 
if precisely chosen, should coincide with the initial values of the physical modes (i.e. the 
Schrodinger modes). Clearly, we would be more happy if these coordinates constitute a 
canonical basis for the reduced phase space. 

We then considered the model of a massive string in a background B-field as an in- 
teresting and nontrivial example for applying the idea of considering boundary conditions 
as Dirac constraints as well as using the symplectic approach for quantizing the model. 
Besides characteristics such as arising mixed boundary conditions (Eqs. [I5Jand[Tn]) and the 
appearance of two sets of infinite number of constraints at each boundary (Eq. IT71) . which 
is common with the massless case, the problem of massive bosonic string has its special 
attractions for two following reasons; 

i) The boundary conditions (115]) and ([TBI) incorporate coordinate and momentum fields 
almost on the same footing. This results to expansions ( |23l) and ( |25l) . in which canonical 
conjugate pairs are present both in X(a, r) and P(<J, r), which finally leads to noncommu- 
tative momentum fields (see Eq. 134]) as opposed to the massless case. 

ii) Upon imposing the constraints, the massive bosonic fields acquire nontrivial zero 
modes (see Eqs. [251) which can not be derived from the limiting case of oscillatory modes. 
We suggest to define the zero mode as solution of generalization of the boundary conditions 
to the whole medium, instead of the boundaries alone. Such a solution clearly satisfies 
the required conditions at the boundaries and should be included in the most possible 
expansions of the fields. We showed that this inclusion could be happen naturally and need 
not to be added by hand (see our discussion after Eq. I2T1) . 

Our main result for this model is finding a canonical basis for the reduced phase space 
with a canonical algebra given in ( 1281) . which apart from being physically meaningful, is 
much simpler to work with, compared to that of reference [ID] . We showed also that the 
Hamiltonian of the system is simply the superposition of harmonic oscillators constructed 
over these coordinate-momentum pairs. 

Giving all technical details, we showed that the coordinate and momentum fields are 
non-commutative at the end-points (Eqs. and |4"TI) ; while they remain, after imposing 

the constraints, canonically conjugate to each other in the bulk of string (Eqs. 1551 and 135]) . 
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